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Abstract: We propose a group theoretic condition which may be applied to extensions 
of the Standard Model in order to locate regions of parameter space in which the elec- 
troweak phase transition is strongly first order, such that electroweak baryogenesis may be 
a viable mechanism for generating the baryon asymmetry of the universe. Specifically, we 
demonstrate that the viable corners of parameter space may be identified by their prox- 
imity to an enhanced discrete symmetry point. At this point, the global symmetry group 
of the theory is extended by a discrete group under which the scalar sector is non-trivially 
charged, and the discrete symmetry is spontaneously broken such that the discrete sym- 
metry relates degenerate electroweak preserving and breaking vacua. This idea is used to 
investigate several specific models of the electroweak symmetry breaking sector. The phase 
transitions identified through this method suggest implications for other relics such as dark 
matter and gravitational waves. 
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1 Introduction 

Standard cosmology of the early universe within the context of a large class of models 
embedding the Standard Model (SM) of particle physics predicts the existence of an elec- 
troweak symmetry breaking (EWSB) phase transition (PT). Collider constraints alone 
cannot determine the nature of the EWSB PT in a model independent way. However, 
additional information is available in the form of cosmological relics, which were produced 
in the early universe and survive as direct probes of the physics of the era during which the 
temperature was electroweak scale. Relics such as the baryon asymmetry [1], primordial 
gravitational waves [2—5], and (modifications to) the dark matter relic abundance [6-9], 
may have been generated at the electroweak scale PT(s). 

Generating the baryon asymmetry through CP violations at electroweak symmetry 
breaking bubbles [1], requires a strongly first order phase transition (SFOPT) to protect 
the baryon number in the broken phase. In this context, a SFOPT may be defined as a 
first order PT in which the (thermal) expectation value of the SM-like Higgs v(T) = (h) 
satisfies v(T)/T > 1 in the broken phase after the phase transition completes, such that 
weak sphaleron processes are inactive [1, 10]. It is well-known that the SM is unable 
to accommodate a SFOPT while satisfying the Large Electron-Positron (LEP) Collider 
bounds on the Higgs mass [11]. This is one of the main motivations for considering an 
extended Higgs sector. Many beyond the Standard Model theories are able to accommodate 
a SFOPT, including supersymmetry, two Higgs doublet models, and minimal scalar singlet 
extensions of the SM. However, if the extra scalar fields obtain vacuum expectation values 
(vevs), one often finds that new patterns of symmetry breaking become accessible. This 
fact makes the phase transition more difficult to study, because quantities such as v(T)/T 
are nonanalytic functions of the parameters of the model. Consequently, many beyond the 
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Standard Model PT analyses rely on an intensive numerical parameter scan to search for 
SFOPT. Although such scans may be capable of locating SFOPTs, on their own they do 
not reveal why one particular parametric limit is favored over another. 

In this article, we propose a group theoretic guideline which will aid the search for 
SFOPT in a large parameter space and help to identify why certain parametric limits are 
favored over others. Our guideline is motivated by the following heuristic argument. In 
perturbative thermal effective potential computations, the thermal mass is of the order 
cT 2 where c is a thermal loop factor. Therefore, if all the renormalized coupling constants 
are of order unity and all mass scales are of the electroweak scale, we expect that the phase 
transition will occur at a temperature T ~ vj 'i/c such that v{T)/T ~ y/c < 1, and the PT 
is typically not strongly first order. Hence, in order to have a SFOPT, the renormalized 
parameters of the theory must be near a special point in the parameter space. An ideal 
parametric limit which overcomes the natural thermal loop suppression is the region where 
v(T)/T — > oo. To achieve this, it would be unnatural to expect v(T) to deviate by many 
orders of magnitude from the electroweak scale, because of the constraint that v(0) defines 
the electroweak scale. On the other hand, v{T)/T may be enhanced by taking the T — ¥ 
limit. 

The limit of low phase transition temperature and large v(T)/T can be achieved nat- 
urally by employing a discrete symmetry. The phase transition begins at the critical 
temperature T c , defined as the temperature above which the thermal corrections are suffi- 
ciently large as to make the EW symmetric phase energetically favored, and below which 
the EW broken phase is favored. Hence, at T = T c the thermal effective potential possesses 
two degenerate minima corresponding to the EW symmetric and broken phases (see also 
Appendix A). One may enhance v(T c )/T c by taking T c — > provided that there is a mech- 
anism guaranteeing that the theory possesses such degenerate vacua even in the absence 
of thermal corrections. One mechanism that yields degenerate vacua is the spontaneous 
symmetry breaking of a discrete group (see e.g. [12, 13]). After spontaneous symmetry 
breaking, one finds a set of degenerate vacuum states which fall into a coset representation 
of the discrete group. Moreover, if the discrete symmetry group does not commute with 
the electroweak group, then the scenario described above may be achieved: the electroweak 
symmetry is broken in one vacuum and unbroken in a second degenerate vacuum implying 
T c = and v(T c )/T c = oo. 

Of course the existence of degenerate vacua alone does not imply v(T)/T 3> 1, since 
the EW phase transition must take place, and this is not necessarily the case in extensions 
of the SM with multiple vacua. If the discrete symmetry is exact, then T c = and the phase 
transition does not proceed because the broken phase never becomes energetically favored. 
Hence, we will consider models in which the discrete symmetry is generally approximate, 
but becomes exact at a particular parametric point, referred to as an enhanced discrete 
symmetry point (EDSP) 1 . Then the heuristic arguments above imply that one can expect 
to find SFOPT in a parametric neighborhood of an EDSP and connected to it by a con- 

1 In general, a model may possess multiple EDSPs each relating the EW broken and symmetric vacua 
by a different symmetry transformation. 
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tinuous "small" deformation which breaks the discrete symmetry. Precisely how "small" 
a deformation is required depends upon two model-dependent conditions: the condition 
that the electroweak PT must complete and upon the order unity number that sits at the 
right of the inequality v(T)/T > 1. Hence the takeaway message is that one can make the 
analysis of and search for SFOPT in a large parameter space more tractable with the aid 
of an EDSP "lamppost" which signals the parametric neighborhood which is favorable for 
SFOPT. 

The order of presentation is as follows. In Sec. 2 we motivate our group theoretic 
identification of SFOPTs. In Sec. 3, we employ our technique to explore three example 
models. We then finish with some concluding remarks in Sec. 4 and an appendix which 
reviews some relevant basics of phase transitions used in this paper. 

2 Why Discrete Symmetry? 

Suppose that a given theory is exactly invariant under an internal discrete symmetry group 
G. It is well-known that the spontaneous symmetry breaking of G down to H C G leads 
to the vacua giving a nontrivial coset G/H representation [12, 13]. We will first illustrate 
how this connects to a SFOPT in a perturbative single real scalar field toy model, and then 
proceed to give a more general discussion. 

Quite often in extensions of the SM, other scalar fields along with the Higgs obtain vevs 
at the electroweak phase transition. One may model such a first order phase transition with 
the following toy theory in which ip represents the linear combination of the SM Higgs and 
other scalar fields. Consider the theory of a real scalar field <p with the classical potential 



and suppose that cp is coupled to a family of N fermions C D {mi + hnp)ipiipi. Note 
that this theory has no internal symmetries for non-special values of the parameters 
{M 2 ,£, X,hi,mi}. When we turn on temperature, there will be a thermal bath of <p and 
tpi particles. If the fermions are relativistic at the electroweak scale (i.e., mf <C T 2 ), then 
the thermal effective potential can be written to leading order as 



neglected the subdominant thermal corrections (such as the non-analytic term) and the h 
radiative corrections. 

As long as the supercooling is small (e.g., as measured by the fractional temperature 
change during the duration of the PT), the PT occurs at the temperature near T c at which 
the thermal effective potential V e s displays two degenerate minima (for more details, see 
Sec. A). Solving this constraint for T c gives 



(2.1) 



y eff ( v ,T)^%) + c rV 



(2.2) 



where c ~ Ntif/12 [14]. Here, in the so-called high-temperature approximation, we have 




(2.3) 
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In this simple toy model and subject to the approximation Eq. (2.2), there is an enhanced 
Z2 symmetry at T = T c . Explicitly, the potential at T = T c becomes 



\<p 2 
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V eS (ip,T c ) 



4 



(2.4) 



and respects the discrete symmetry 



( 



-Sh-K 



) 



(2.5) 



which was not originally present in Eq. (2.1). This Z2 exchanges the degenerate vacua at 
<p = and ip = v(T c ) = 2£/X across a potential barrier at 92 = £ /A. It is to be noted that 
Eq. (2.4) is independent of M, and thus this symmetry exists in this toy model for any 
critical temperature T c that can be tuned using M? 

Although there is no electroweak symmetry in this toy model, there is still a first order 
phase transition, and we can investigate the parametric dependence of its order parameter 
v(T c )/T c . Since v(T c ) = 2£/Xis independent of M 2 , the order parameter can be maximized 
by varying M 2 to minimize T c . Even though the high-temperature expansion breaks down 
when T drops below the mass of the fermion, the formal limit T c — > can be taken assuming 
that the fermions are massless. The formal solution to T c = is 3 



The important observation is that 1—a = corresponds to an EDSP in the parameter space 
at which the zero-temperature scalar potential Eq. (2.1) is invariant under the symmetry 
transformation Eq. (2.5). At the EDSP, 1 — a = 0, the order parameter 



formally diverges, and for 1 — a <C 1, the phase transition may be made arbitrarily strongly 
first order. 

Hence, our group theoretic guideline leads us to identify the parametric region in the 
vicinity of the EDSP 1 — a = as favorable for SFOPT. However, for this region to be truly 
viable, it must be the case that the rate at which bubbles of the broken phase nucleate is 
sufficiently large that the phase transition actually completes. This requires the discrete 
symmetry to be weakly broken, such that the PT occurs at a nonzero T. In the toy 
model, such breaking can be accomplished explicitly at the classical level through a finite 
excursion from the EDSP (i.e. 1 — a = e 7^ 0), or radiatively through the Yukawa coupling. 
Indeed, in many extensions of the SM where singlets are introduced, the relevant discrete 
symmetry transforms both the Higgs and the singlet fields. Since the singlets lack SM 

2 Although this is an enhanced Z2 symmetry at T = T c , since the symmetry does not generically exist 
at other temperatures, it is not the enhanced discrete symmetry point relevant for this paper. See below 
for further clarification. 

3 The idea of focusing on T c ~ was recently emphasized by [15]. 

4 Note that the bubble nucleation rate is zero at T = T c . 



a = \M 2 /2£ 2 = 1 . 
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Figure 1. The order parameter, calculated as v(T c )/T c (dashed) and v(T r )/T r (solid), plotted 
against a = \Al 2 /2£ 2 . The insets show U(tp) for particular values of a in each of the associated 
colored regions. 



gauge couplings, radiative corrections necessarily break the discrete symmetry to a degree 
controlled by the strength of the gauge interactions. If the breaking of the symmetry is 
so large that the potential does not have the qualitative features of Eq. (2.1) near a = 1, 
then the EDSP method loses its advantage for identifying SFOPT. If the breaking of the 
symmetry is so small that bubbles will not nucleate fast enough to complete the PT, then 
any candidate parameter points found with the EDSP method are inherently not viable. 
Since this non-completion of the PT will be a general feature of the region of parameter 
space nearby to the EDSP, we must take extra care in choosing the temperature at which 
to evaluate the EW order parameter v(T)/T. Up to this point in the discussion, we have 
evaluated v(T)/T at the degeneracy temperature T c , which, physically, is the temperature 
in the symmetric phase at the onset of supercooling. However, a first order phase transition 
proceeds with the nucleation of bubbles of broken phase which subsequently collide and 
reheat the plasma to a temperature T r (see App. A for precise definition). Since the purpose 
of the SFOPT criterion v(T)/T > 1 is to ensure suppression of weak sphaleron processes 
in the broken phase after the phase transition, the most physically relevant temperature at 
which to evaluate v(T)/T is the reheat temperature T r . 

To obtain a numerical intuition for our proposal, consider Fig. 1 where we have plotted 
v(T c )/T c (dashed) and v(T r )/T r (solid) while varying a and fixing U'(v) = at v = 
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300 GeV, U"(v) = (50 GeV) 2 , N = 1, and h = 0.3. In this figure, we also show U(<p), such 
as to make the discrete symmetry evident at the EDSP. As expected, v(T c )/T c diverges at 
the EDSP and is arbitrarily large for arbitrarily small discrete symmetry breaking (1 — a <C 
1) . On the other hand, v(T r )/T r cannot be calculated if the discrete symmetry is too weakly 
broken (1 — a < 0.5), because the phase transition does not occur. However, sufficient 
discrete symmetry breaking (1 — a > 0.5) yields SFOPT which become monotonically 
weaker as the degree of symmetry breaking grows. We have used the same coloring in 
Fig. 1 as we do in the rest of this article to distinguish the varous regions of parameter 
space: the phase transition does not occur because the broken phase is not energetically 
favored (green); the PT does not occur because the bubble nucleation rate is too low 
(orange); a strongly first order PT occurs (blue); a weakly first order or second order PT 
occurs (gray); the EDSP (purple dot); and the point at which the barrier disappears (red 
dot). 

Now let us return to a more broad discussion of the connection between discrete sym- 
metry and strongly first order phase transition. In retrospect, we recognize that the ex- 
istence of an EDSP associated with a discrete symmetry under which the vacua form a 
coset representation (along with the condition that spontaneous symmetry breaking oc- 
curs) is sufficient to obtain v(T c )/T c — > oo since T c = implies a degeneracy at the level 
of the nonthermal effective potential. Even though the toy model calculation was accom- 
plished using the leading high-temperature T dependence and the classical potential, this 
statement regarding the EDSP is an exact statement for an exact effective potential. In 
other words, as far as this exact statement is concerned, it is not particularly important 
that T = T c corresponded to an enhanced symmetry point for general T c as in the case 
of this simple one dimensional toy model (see Eq. (2.4)), nor is it important that quan- 
tum radiative corrections from the Yukawa couplings break the discrete symmetry given 
by Eq. (2.5). One final ingredient, which is important for electroweak baryogenesis but is 
not represented in the toy model is that at least two vacua in the coset space must carry 
different electroweak quantum numbers. Otherwise, the PT will not be an electroweak 
symmetry breaking PT. This means that the discrete group must not commute with the 
electroweak group and one element in the coset representation must be an electroweak 
singlet. Hence our group theoretic guideline may be summarized as: an arbitrarily strong 
phase transition (i.e., v(T c )/T c ^> 1) may be found in the parametric neighborhood of an 
EDSP if 1) the condition for spontaneous (discrete) symmetry breaking is satisfied (such 
that there will be degenerate vacua), 2) the discrete group does not commute with the 
electroweak group, and 3) its coset representation contains an electroweak singlet element 
(such that the EW symmetry is broken in one vacuum and preserved in another). 
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3 A Few Examples 
3.1 SM with Low Cutoff 

As a first example, we will consider a generic extension of the SM with a low scale cutoff, 
as studied by [16—18]. Provided that the UV physics does not violate the EW symmetry, 
then upon integrating it out one obtains a classical potential of the form 

- £ 3A(|^|4) 2 + i,(|«t ff |_^ 3 (3.!) 
up to terms of order H 8 /A 4 . Writing the Higgs doublet in terms of the fundamental scalar 

T 

Higgs h as H = (0, h/y/2) , and using m 2 H = 2Xv 2 , the potential becomes 

up to constant and higher order terms. There exists an enhanced discrete symmetry point 5 
at which a Z 2 symmetry is nonlinearly realized, 

EDSP: m H A = v 2 Z' 2 : h ->• ~ + J v 2 - ^h 2 . (3.3) 

The 7U 2 symmetry exchanges the minima at h = and h = v while leaving the maximum at 
h = u/\/3 invariant. We have reproduced an earlier PT analysis [16] in order to illustrate 
the proximity of SFOPTs to the EDSP. Moreover, we have extended the previous analysis 
by calculating the more physically relevant order parameter v(T r )/T r , instead of v(T c )/T c . 
Our results are summarized in Fig. 2, and are in good agreement with Fig. 2 of [16] which 
shows the same slice of parameter space. We find that nearby to the EDSP (purple curve), 
the PT is strongly first order (blue), and that the PT becomes weaker moving away from 
the EDSP. It is also worth noting that while the barrier persists, the PT most likely does 
not occur, as evidenced by the lack of blue in the region between the purple and red curves 
except for a small sliver above mu = 200 GeV. 



3.2 SM Plus Real Singlet xSM 

Next, we will consider models with multiple scalars in the electroweak sector. Extending 
the SM by a real scalar singlet s, we obtain a model known as the xSM [19], which has the 
classical potential 

U(h, s) = ^ - ^h 2 + ^ 4 + ^ 3 + K 2 + % s 2 h 2 + ^sh 2 . (3.4) 
4 2 432 4 4 

5 It may be more appropriate to use the term "enhanced discrete symmetry plane," as the condition 
mfrA = v 2 actually specifies a hypersurface in the parameter space, but we will continue using EDSP for 
simplicity. 
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Figure 2. The parameter space nearby to the EDSP (purple curve). The coloring is the same 
as in Fig. 1. The PT order parameter v(T r )/T r is indicated by the overlaid contours. SFOPTs are 
found in the blue region and become weaker in the gray region, farther from the EDSP. 

Since there is no symmetry protecting s = 0, generally both h and s will obtain vevs, 
denoted v and xq respectively, and the mass parameters may be written as 

/j 2 = X oV 2 + °^-xl + ^-x and b 2 = -hx% - b 3 x - -^-v 2 - . (3.5) 

II I 4 Xq 

Provided that x$ ^ 0, the cubic terms s 3 and sh 2 help to generate a barrier separating 
the symmetric and broken vacua and make the PT strongly first order. A number of PT 
analyses [15, 20-22] have revealed that the xSM can accommodate a strongly first order 
electroweak PT. They also find that this model displays multiple patterns of symmetry 
breaking such that, either h and s can obtain vevs at the same temperature, or s can receive 
a vev prior to electroweak symmetry breaking. If we were to search for SFOPT by randomly 
choosing order one parameters, there would be no way of anticipating what pattern of 
symmetry breaking would be realized, or if the EW symmetry would be spontaneously 
broken at all. Moreover, since Eq. (3.4) has six free parameters, such a random search 
could become quite time consuming. 

The discrete symmetry technique greatly simplifies the SFOPT search. We are able to 
specify a desired pattern of symmetry breaking to investigate, identify the corresponding 
discrete symmetry, compute the associated EDSP, and begin searching by perturbing from 
the EDSP. Here, we will focus on a particular pattern of symmetry breaking in which both 
s and h obtain vevs simultaneously, and we will compare our calculation against the "high- 
T trivial singlet vev" case of [20]. The appropriate discrete symmetry is a Z2 relating the 
vacua at {h, s} = {0, 0} and {v, xo}. We can identify the associated EDSP by first reducing 
Eq. (3.4) to Eq. (2.1) and then imposing a = 1. This is accomplished by focusing on the 
one-dimensional linear trajectory {h, s} = {v, xq} x <p/yj v 2 + Xq parametrized by ip, which 
interpolates between the EW-symmetric and EW-broken vacua. Along this trajectory, the 
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Figure 3. A slice of the xSM parameter space showing the proximity of SFOPT (blue region) to 
the enhanced symmetry axis (purple curve). 

potential can be written in the form of Eq. (2.1) with 

A= \ oV ± + b 4 x* + a 2 v 2 x 2 
{v 2 + xl) 2 

M 2 = ^v 2 + x 2 (?>£ - \^v 2 + x 2 

Then, upon resolving the condition a = 1 we find the enhanced discrete symmetry point, 
EDSP : = 12a 2 v 2 xl + 3aiv 2 x + 4& 3 x[j + 12b±x% + 6A v 4 

: („-§)-»-(„-§). (3.7) 

In general, the PT will not occur along the trajectory parametrized by <p, but nevertheless 
this linear interpolation is useful for identifying the EDSP. 

Once again, we have numerically investigated the strength of first order PTs in the 
vicinity of the EDSP. We have chosen a parameter set which allows us to reproduce Fig. 4 
(left panel) of [ ] by fixing a x = -933 GeV, a 2 = 0.69, 6 3 = 356 GeV, b 4 = 0.53 and 
scanning Ao £ [0,1] and log 10 xo G [—1,3]. Our results are shown in Fig. 3. A few 
observations may be made. First, as anticipated, the first order PTs are strongest close to 
the EDSP (purple) curve and become weaker farther away. Second, there is a large region 
(green) in which the EW remains unbroken. Below the EDSP (purple) curve, the origin 



a 



x (3aiu 2 + 463x5) 
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(3-6) 
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remains the global minimum of the effective potential, whereas at large values of xq > 10 2,1 , 
the global minimum sits at s < 0. Third, in comparing with [20], one must bear in mind 
that we have fixed the remaining parameters, whereas those authors have scanned the full 
parameter space and projected onto these coordinates. As such, the region where we find 
SFOPT is much smaller than what is suggested by Fig. 4 of [20]. However, this just goes to 
show that it is typically difficult to find SFOPT in a large parameter space without either 
a large parameter scan or some guiding principle. 



3.3 SM Plus Real Z 2 -Charged Singlet Z 2 xSM 

As a final example, we turn out attention to the Z2XSM, which extends the SM by a real 
scalar singlet s such that the scalar potential becomes [15, 23] 

U(h, s) = -hr - — hr + — s* + — s z + — s'h- . (3.8) 
4 2 4 2 4 

The singlet is charged under a Z2, which restricts the allowed operators, but extends the 
possible patterns of symmetry breaking, because now (s) = is radiatively stable. We 
will focus on a particular parameter region in which there is transitional Z2 symmetry 
breaking: at temperature T > T a both Z2 and the EW symmetry are restored, at T = T a 
the singlet obtains a vev breaking Z2, and at T = TJ, < T a the Higgs field obtains a vev 
and the singlet's vev returns to zero, thereby breaking the EW symmetry and restoring the 
Z2 (i.e., EW x Z 2 — > EW x j^f — >£Wx Z2). In the context of this pattern of symmetry 
breaking, the enhanced discrete symmetry point admits an §> 2 symmetry, 

EDSP : 64 = A and b 2 = -fi 2 §2 : h s (3.9) 

where we will also take 02 > 2A to ensure that the discrete symmetry interchanges vacua. 
Note that this §2 symmetry is more restrictive than the Z2 symmetries we considered in 
the previous examples. To illuminate the role of the EDSP in locating SFOPT, we will 
reparametrize 64 = A + A64 and 62 = —/J- 2 + A62 to write the potential as 



uu,*)= ^(/ i 4 + s 4 )-^(/ i 2 + , 2 ) + |/,V 



+ 



A64 4 A6 2 2 

— H — —s 



(3.10) 



where we have also used /i 2 = At> 2 . In this parameterization, we expect to find SFOPT 
nearby to the EDSP at A6 4 = A6 2 = 0. 

We present the results of our numerical analysis in Fig. 4, where we have fixed A ~ 0.12 
to give a Higgs mass 6 of rrih = V2Xv 2 = 120 GeV. As in the previous examples, the phase 
transition strength decreases monotonically with distance from the enhanced symmetry 
axis. Significantly far from the EDSP, the phase transition proceeds with a different pattern 
of symmetry breaking. In the brown region, the EW symmetry breaks without transitional 



6 Since the axes of Fig. 4 depend only on the ratios A64/A and A&2/A?; 2 , a change in the Higgs mass (via 
A) could be absorbed by A64 and A&2, such that the qualitative features of Fig. 4 would remain unchanged. 
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Figure 4. Three slices of the Z2XSM parameter space for fixed A 0.12. The origin A62 = A64 = 
is an EDSP at which the theory has an §2 discrete symmetry. 

Z2 violation (EW x Z2 — t--EWx Z2), in the yellow region the Z2 remains broken in the low 
temperature vacuum (EW x Z2 — > EW x^f— t--EWx and in the purple region there 
exists an intermediate phase in which both Z2 and the electroweak symmetry are broken 
(EW xZ 2 4 EWxX^^xX^^x Z 2 ). 

The region of parameter space nearby to the EDSP displays an interesting phenomenol- 
ogy. Since the singlet mass is given by 



one typically finds m s < rrih nearby to the enhanced symmetry point. The unbroken 
Z2 symmetry ensures that the singlet is stable, and thus it is a dark matter candidate 
which annihilates to Higgses with a cross section proportional to o|. A number of analyses 
[19, 24-31] have considered this scenario and found that a 2 and the singlet mass m s can 
be strongly constrained by assuming that the s particle composes all of the dark matter. 
Collider experiments, such as the LHC, may also be able to constrain the Higgs-singlet 
coupling. For A6 2 /A-u 2 < (3 — a2/A)/2, the singlet mass is less than half of the Higgs 
mass and the invisible decay channel h — > ss becomes kinematically accessible. Then, 
a measurement of the invisible decay width may constrain the Higgs-singlet coupling 02 
[18, 31-34]. On the other hand, the singlet self-coupling 64 remains unconstrained. This 
is because unlike in other limits of this model and similar models [32, 33, 35, 36], the 
unbroken Z2 symmetry prevents the Higgs and singlet from mixing. Consequently, the 
singlet self-coupling 64 is practically impossible to constrain at colliders, and contributions 
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to the anomalous Higgs trilinear coupling [37] are loop suppressed. Finally, let us point 
out that the transitional Z2 violation limit may not suffer from the domain wall problem 
that generally accompanies models with spontaneously broken discrete symmetries. When 
the Z2 breaks in the first step of the PT, domain walls will be generated. However, once 
the EW symmetry is broken and the Z2 symmetry is restored, the domain walls should be 
"wiped out" by the Z2-symmetric vacuum field configuration. This may lead to a unique 
gravitational wave spectrum. 

4 Conclusion 

Strongly first order phase transitions (SFOPTs) are required for electroweak baryogenesis 
and may have other interesting implications for early universe relics. In this article we 
have discussed a general analytic guideline, based on symmetry principles, which is useful 
in identifying a region of parameter space favorable for SFOPT: an arbitrarily strong PT 
can be found for parameters near an enhanced discrete symmetry point (EDSP) if the 
condition for spontaneous symmetry breaking is met and if the discrete symmetry relates 
the electroweak symmetry preserving vacuum to one in which it is broken. Group theoreti- 
cally, this means that the coset representation of the broken discrete symmetry contains an 
electroweak singlet and the discrete group does not commute with the electroweak group. 
Because of phenomenological requirement of completing the PT at a nonzero tempera- 
ture, the symmetry must be broken by parametric deformations away from the EDSP. 
As the deformation decreases, the strength of the PT tends to increase. We applied this 
guideline to study the electroweak PT in three specific models. In each of the models 
considered, SFOPTs occur in close proximity to the EDSP, as expected. In this way, the 
enhanced symmetry point acts like a lamppost in the parameter space, signaling the loca- 
tion of SFOPTs. It would be interesting to apply a similar EDSP-motivated analysis of the 
electroweak phase transition to models with larger scalar sectors and greater parametric 
freedom, such as singlet extensions of the Minimal Supersymmetric Standard Model. 

It is not unnatural to expect SFOPT to be localized in the vicinity of an EDSP. Strongly 
first order phase transitions almost always require some fine-tuning of the parameters in 
the theory. From an UV completion point of view, such fine-tuning could be more natural 
if it is close to a point of the parameter space with enhanced symmetry. It is also clear that 
degenerate vacua may be found even without discrete symmetry, and thus our guideline 
provides a sufficient, though not necessary, condition for locating SFOPT. Nonetheless, 
such parametric regions form a large class of possibilities which can most likely always 
occur in practice. 

We also observe (as did [7, 15]) that the PT tends not to proceed at all unless the barrier 
separating the EW-broken and EW-unbroken vacua is very small or not present at all (along 
the red curve), because otherwise the tunneling rate is too strongly suppressed. Hence, 
the deformations away from the EDSP required for phenomenologically viable SFOPTs 
are not vanishingly small and are model dependent. Although such phenomenologically 
viable parametric regions can be arrived at by deforming away from enhanced continuous 
symmetry points rather than EDSPs, the EDSP starting point guarantees the existence 
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of potential barriers required for a first order PT. In that sense, our proposal here is 
advantageous over the enhanced continuous symmetry point perspective. 

Proximity to an EDSP implies interesting relations between parameters in the extended 
Higgs sector, which is responsible for the dynamics of the electroweak symmetry breaking. 
Such relations will manifest themselves in both the spectrum of the states in the Higgs sector 
and their couplings. Probing this sector is the central scientific focus of the LHC. We might 
have already seen the discovery of the Higgs boson on the horizon [38, 39]. Discovering 
the additional states in the extended Higgs sector and measuring the parameters in the 
Higgs potential are expected to be very challenging tasks. At the same time, confirming 
the structure of the Higgs sector to be consistent with a SFOPT would establish a striking 
link to the generation of the baryonic asymmetry in the universe. 

A Details of Phase Transition Calculation 

For the phase transition analyses in this paper, we have calculated the thermal effective 
potential V e fi(<p,T) through one-loop order using the standard techniques [40-42]. We 
numerically minimize V e g with respect to <j) to obtain the scalar field expectation values in 
the symmetric and broken phases, v sym (T) and Vbrk{T), respectively. The latter quantity 
is sometimes referred to in the text as simply v(T). The critical temperature T c is defined 
as 

V eS (v sym {T c ),T c ) = V cS (v hrk (T c ),T c ) . (A.l) 

We use V e s(cj),T) to calculate the action 8 Ss(T) of the bubble field configuration that me- 
diates the vacuum transition [45-48] . We determine the bubble nucleation temperature T n 
by requiring the bubble nucleation rate per Hubble volume to exceed the Hubble expansion 
rate. This condition may be resolved to 

S 3 (T n )/T n = 140 (A.2) 

where the value on the right hand side depends only logarithmically on the model pa- 
rameters [49, 50]. Finally, we calculate the temperature T r of the plasma after the phase 
transition ends and the plasma has been reheated. This is obtained by assuming that the 
universe does not expand significantly during the phase transition and then by imposing 
energy conservation [ ] 

Psym(T n ) = Pbrk(^r) (A.3) 

where 

p(T) = V cS (v(T),T) - T-^V eS (v(T),T) (A.4) 

is the energy density in the symmetric or broken phase, respectively. 

7 This definition of v(T) implies that T c ,T n , and T r will be dependent upon the choice of gauge [43, 44]. 
Though this may affect the numerical accuracy of our results, we expect that the qualitative parametric 
dependence of the EW order parameter nearby to an EDSP, which is our primary interest, will remain 
unchanged. 

8 For the models of Sees. 3.2 and 3.3 which have more than one scalar field participating the phase 
transition, we calculate the bounce using the approximation described in [6]. 
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